Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 1 February 2008 (MN MfeK style file vl.4) 



X-ray irradiation in low mass binary systems 



Guillaume Dubus 1 , Jean-Pierre Lasota 1 , Jean-Marie Hameury 2 and Phil Charles 4 

1 UPR 176 du CNRS, Departement d'Astrophysique Relativiste et de Cosmologie, Observatoire de Paris, Section de Meudon, 
F-92195 Meudon Cedex, France, dubus@obspm.fr, lasota@obspm.fr 

2 UMR 7550 du CNRS, Observatoire de Strasbourg, 11 rue de I'Universite, F-67000 Strasbourg, France, hameury@astro.u-strasbg.fr 
4 Astrophysics, University of Oxford, Nuclear and Astrophysics Laboratory, Keble Road, Oxford OX1 3RH, UK, pac@astro.ox.ac.uk 



OO 
On 

o\ 

Or 
m 



> 

m 
o 

o 

oo 

On 



1 February 2008 



ABSTRACT 

We calculate self-consistent models of X-ray irradiated accretion discs in close binary 
systems. We show that a point X-ray source powered by accretion and located in the 
disc plane cannot modify the disc structure, mainly because of the self-screening by 
the disc of its outer regions. Since observations show that the emission of the outer 
disc regions in low mass X-ray binaries is dominated by the reprocessed X-ray flux, 
accretion discs in these systems must be either warped or irradiated by a source above 
the disc plane, or both. We analyse the thermal- viscous stability of irradiated accretion 
discs and derive the stability criteria of such systems. We find that, contrary to the 
usual assumptions, the critical accretion rate below which a disc is unstable is rather 
uncertain since the correct formula describing irradiation is not well known. 
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1 INTRODUCTION 

There is a convincing body of evidence that most of the 
optical light emitted by accretion discs in Low Mass X-ray 
Binaries (LMXBs) i s due to reprocessed X-rays illum inat- 
ing the di sc surf ace ( van Paradijs fc McClintock, 1995 ). van 
Paradijs (1996) pointed out that this irradiation, observed 
both in steady LMXBs and during the outburst phases of 
LMXB transient sources ("Soft X-ray Transients", SXTs) 
should influence the stability properties of accretion discs in 
such systems. X-ray irradiation increases the surface (and 
therefore the central) disc temperature and could, in prin- 
ciple, prevent the growth of the thermal instability that ap- 
pears at a temperature for which hydrogen becomes partially 
ionized. This instability leads, under certain conditions, to 
a global thermal- viscous instability which is taken to be re- 
sponsible for dwarf nova outbursts and SXT events (see e.g. 
Cannizzo 1993). 

Models used to describe X-ray irradiated discs are of- 
ten, as we show below, inconsistent and/or assume unveri- 
fied hypotheses. Despite these drawbacks, results obtained 
in such studies seem to be in agreement with observations. 
In the present work, we construct self-consistent models of 
irradiated disc vertical structures and show what should be 
modified in the usual treatment to account for observations 
whilst being consistent. 

This article is organized as follows: in Section 2 we re- 
discuss the basics of irradiated disc models. In Section 3 we 



present self-consistent models of irradiated disks and discuss 
their stability. 



2 X-RAY IRRADIATION OF ACCRETION 
DISCS IN CLOSE BINARY SYSTEMS 

2.1 Steady discs 

In most studies of X-ray illuminated accretion discs in close 
binary systems it is assumed (explicitly or implicitly) that 
the disc is isothermal in regions in which reprocessed light 
is an important part of the disc emission (see e.g. Vrtilek et 
al., 19 90). For example, de Jong, van Paradijs & Augusteijn 



( |l996| ) use the Vrtilek et al. (1990) model to deduce var- 
ious properties of X-ray irradiated discs in LMXBs, such 
as the albedo and the disc opening angle. It is sometimes 
said that the values obtained in this way are suggested by 
observations, and they are used even i n the case were it i s 
stated that the disc is not isothermal (King et al., 1997a). 
To justify the assumption of an isothermal disc it is argued 
that if the irradiating flux is larger than the flux due to vis- 
cous dissipation in the disc, the viscous heating term can be 
neglected in the energy equation so that the central temper- 
ature is equal to the irradiation temperature. This reasoning 
is, in general, erroneous: irradiation dominates the internal 
disc structure only if the irradiation flux divided by the disc 
total opti cal depth is larger than the flux due to viscous dis- 
sipation ( Ly^Udj^Snrryaev^JljJTfi |_ rucliman, Mineshige & 



Wheeler, 1990t fHubeny, 199lj [Hure et al., 1994| ). This is a 
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simple consequence of energy conservation and of the radia- 
tive transfer equation and is independent of the details of 
disc models. Since there seems to be some confusion about 
this point, we shall derive the relation between disc midplane 
and surface temperatures from the basic equations. 

Since we are interested here in stable and therefore hot 
accretion discs, the convective energy transport does not 
play an important role (even though it is taken into account 
in the numerical model described in section 3.1). We write 
the energy conservation as : 



nj— = Qvis{R, z) 
dz 



(1) 



where F is the vertical (in the z direction) radiative flux and 
Q v is(R, z) is the viscous heating rate per unit volume (R is 
the radius). Eq. (Q) states that an accretion disc is not in 
radiative equilibrium, contrary to a stellar atmosphere. For 
this equation to be solved, the function Q vlB {R,z) must be 
known. Of course, despite many efforts and recent progress 
in understanding the origin of viscosity in accretion discs 
(see e.g. Zahn 1991; Balbus & Hawley 1998) the func- 
tion Q v i s (R,z) is largely unknown. The ansatz most widely 
used in practical applications is based on th e "a - viscos- 
ity" prescription (Shakura & Sunyaev 1973) according to 



which the kinematic viscosity coefficient can be written as 
v — (2/3)aCg /Qk, where a is the viscosity parameter (< 1), 
f^K is the Keplerian angular frequency and c s = \/ P/p is 
the sound speed, p the density, and P the total pressure. The 
viscous dissipation can thereafter be written as (e.g. Smak 
1984) 



Q vls (R,z) = (3/2)afi K P 



(2) 



Viscous heating of this form has important implications for 
the structure of optically thin layers of accretion discs and 
may lead to the creation of coronae and winds (Shaviv & 
Wehrse 1986; 1991). Here, however, we are interested in the 
effects of irradiation on the inner structure of an optically 
thick disc, so the precise form of the viscous heating is not 
important, as ca n be seen by comparing our results with 
those of Hubeny fll99l| ). 

When integrated over z, the rhs of Eq. (|l]) using Eq. (2) 
is equal to viscous dissipation per unit surface: 



Pdz, 



(3) 



which is close, but not exactly equal, to the surface heat- 
ing term (9/8) u'EQ.^ generally used in the literature. The 
difference between the t wo expressions may b e important 
in numerical calculation (Hameury et al., 1998) but in the 
present context is of no importance. 
One can rewrite Eq. (hi) as 

*f = -/(r)^S (4) 
dr Ttot 

where we introduced a new variable, the optical depth dr = 
—KB.pdz, kk being the Rosseland mean opacity and Ttot = 
J^ 00 KRpdz is the total optical depth, /(r) is given by: 



fir) 



(r 



unpdz 



Jo 



Pdz 



KRp 



(•») 



of two rather well defined scale heights, the pressure and the 
opacity scale heights, which are comparable, so that / is of 
order of unity. 

At the disc midplane, by symmetry, the flux must van- 
ish: .F(Ttot) = 0, whereas at the surface, (r = 0) 

F(0) = aT c % = F vis (6) 

Equation (^) states that the total flux at the surface is equal 
to the energy dissipated by viscosity (per unit time and unit 
surface) . The solution of Eq. (^]) is thus 



F(t) = F viB 1 



Ttot 



(7) 



where /J"'"' /(r)dr = Ttot- Given that / is of order of unity, 
putting /(t) = 1 is a reasonable approximation. The precise 
form of /(t) is more complex, and is given by the functional 
dependence of the opacities on density and temperature; it is 
of no importance in this illustrative section which gives the 
physical basis for understanding the exact numerical results 
that will be given in the following sections. We thus take: 

F(r) = F vis (l - JL) (8) 

V Ttot / 

To obtain the temperature stratification one has to solve 
the transfer equation. Here we use the diffusion approxima- 
tion 



(9) 



appropriate for the optically thick discs we are dealing with. 
The integration of Eq. ([]) is straightforward and gives : 



T 4 (t)-T 4 (0) = |t(i 



2Tt< 



Toff 



(10) 



(11) 



The upper (surface) boundary condition is: 

r 4 (0) = iTi+T* r 

where T£ x is the irradiation temperature, which depends on 
R, the albedo, the height at which the energy is deposited 
and on the shape of the disc (see Eq. |l|). In Eq. (jyj) T(0) 
corresponds to the emergent flux and, as mentioned above, 
T c ff corresponds to the total flux (see e.g. Hubeny 1991) 
which explains the factor 1/2 in Eq (|ll[). The temperature 
stratification is thus : 

TV) = ^ ff [T(l-^) + |] + ^ (12) 

which is very similar to the formula obtained by Hubeny 
(1991) for a different viscosity prescription. For Ttot 3> 1 the 
first term on the rhs has the form familiar from the stellar 
atmosphere models in the Eddington approximation. 

For Ttot S> 1, the temperature at the disc midplane is 

3 



Tc = T 4 (T to t) = -TtotTcff + T 4 r 



(13) 



It is clear, therefore, that for the disc inner structure to be 
dominated by irradiation and the disc to be isothermal one 
must have 

Firr = OTL 
Ttot Ttot 



> F vis 



(14) 



As p decreases approximately exponentially, /(t) is the ratio 



and not just F rr 3> F v i s as is usually assumed. The differ- 
ence between the two criteria is important in LMXBs since, 
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for parameters of interest, r to t il 10 2 — 10 3 in the outer disc 
regions . In Section 3.1 we use self-consistent detailed ver- 



tical disc structure calculations in order to determine the 



influence of irradiation on the inner disc structure. Note 
that we neglect here the possible presence of an X-ray ir- 
radiatio n gen erated corona and wind, described by Idan & 
Shaviv (p96f). 



2.2 Irradiation temperature 

Condition <Q is not however a necessary condition for the 
disc emission properties to be dominated by irradiation. It is 
enough that Fin > F v - is for the disc observed luminosity to 
be dominated by reprocessed emiss ion. Observations suggest 
(van Paradijs & McClintock, 1995) that the absolute visual 



magnitude My of nova-like cataclysmic variables which are 



accreting at rates M ~ 2 — 5 x 10 17 g s 1 ( Warner, 1995[) are 



roughly two magnitudes fainter than the Mv's of LMXBs of 
similar size (since one has to compare luminosities and not 
fluxes one should compare systems with the same emitting 
areas). This would imply Fi rl ~ 6-FVi s so that, although the 
optical emission is dominated by reprocessed X-ray radia- 
tion, the inner disc structure is not dominated by irradiation 
(Eq. 13). 

For a point source, the irradiatio n temperature can be 
written as (Shakura & Sunyaev 1973) 



4 _ vMc 2 (l 

irr 



£) -ffirr / d In Hi 



AttoR 2 



R 



V dlnR J 



(15) 



where rj is the efficiency of converting accretion power into 
X-rays, e is the X-ray albedo and H[ lr is the local height 
at which irradiation energy is deposited, or the height of 
the disc "as seen by X-rays". We use here Hi IT and not H, 
the local pressure scale-height, as is usually written in the 
literature because, in gener al, jj~j rr 7^ H. 



In the usual approach (Vrtilek et al. 199C) one assumes 
that the disc is isothermal, which, for a Keplerian disc, 
implies tha t the pressure scale-height varies like i? 9//r . de 
Jong et al. (1996) applied an isothermal model for X-ray re- 
processing in neutron-star LMXBs and concluded that the 
disc opening angle is ~ 12°, i.e. H ll:l (Rd)/Rd ~ 0.2, where 
^irr(-Rd) is measured at the outer disc radius _Rd- It is easy 
to see that the value 0.2 cannot be the pressure scale-height 
to radius ratio since this would imply temperatures > 10 7 
K at the outer disc radius. In addition, the value H/R = 0.2 
was obtained by modeling two dipping sources in which the 
accreting object is a neutron star. It is not at all o bvious that 
this value applies to all LMXBs. de Jong et al. (1996) also 



deduced a high value for the X-ray albedo e 0.9, although 
they discussed possible lower values for this parameter. 

Equation (113) has been extensively used in recent work 
on disc stability properties in LMXBs. van Paradijs (1996) 
suggested that X-ray illumination might modify the condi- 
tions under which accretion disks undergo dwarf nova e type 
thermal- viscous instabilities. King and co l laborators, ( King 



Kolb fe Burderi, 1996; King & Kolb, 1997; King et al, 1997e 



King et al., 1997b) have considered different implications of 
this idea on the evolution of various classes of LMXBs. All 
these models used Eq. Jlj| ) with r\ £s 0.1, e w 0.9. Since the 
stability criterion applies to the outer disc radius the value 
0.2 is a ssume d for Hi rl /R and din _Hj rr /dln R — 9/7. In King 
et al. ( |l997Ej ), the rhs of Eq. (|l|) is multiplied by H/R for 



black hole LMXBs, supposedly to account for the n on-point 
source character of the X-ray emitter in this case (Shakura 



Sunyaev 1973) 



( [19961 ) 



The new stability criterion proposed by van Paradijs 
seems to correspond quite well to observed proper- 



ties of LMXBs in the sense that most sources that should 
be stable according to this criterion are more or less steady 
and all transient sources are unstable. This new criterion is a 
sufficient condition for disc stability (or necessary for insta- 
bility) . The problem is that it is based on a model which can- 
not represent the discs to which it is applied. The isothermal 
disc model gives an incorrect description of X-ray irradiated 
discs in LMXBs. It is often argued that this is just a sim- 
ple disc model. This is true. Unfortunately, like many other 
simple models, it does not apply to X-ray irradiated discs 
in LMXBs. It assumes, contrary to observed properties of 
accretion discs in LMXBs, that the irradiation temperature 
is equal to the midplane disc temperature. The midplane 
temperature at which a hot, steady, accretion disk becomes 
unstable is (this corresponds to the minimum surface density 
for hot stable discs, E m in) 



21700 K a" 



Mi 

M© 



( R ) 

VlO 10 cm J 



(16) 



which for standard parameters (a ~ 0.1) is close to 30 000 
K. Eq. (|l^) is based on numerical calculations by Hameury 
et al. (1998, hereafter H98); similar formulae were obtained 
by other authors (e.g. Cannizzo 1993b) . Isothermal vertical 
structures are obtained only for irradiation temperatures 
corresponding to excessively high accretion rates and/or 
radii. 



2.3 Stability of X-ray irradiated discs 

However, X-ray irradiation does not have to dominate the 
vertical structure to play an important role in t he di sc prop- 
erties. Also, as pointed out by King & Ritter (1998), X-ray 
irradiation may affect the temporal behaviour of transient 
LMXB outbursts. 

It is instructive to see how irradiation modifies disc sta- 
bility properties. The effective temperature at which a hot 
disc becomes thermally unstable is, according to H98, given 
by 



T eff ,B = 7200 a' 



M ) \ 



R 



10 10 cm 



K (17) 



For parameters of interest, this leads to tr, ~ 10 2 (see Eq. 
16) and, as mentioned above, the vertical disc structure can 
then only be changed for very high irradiation fluxes. Much 
lower fluxes, however, are sufficient to modify its thermal 
stability properties, as will be shown in the next sections. 

The disc is thermally stable if radiative cooling varies 
faster with temperature than viscous heating (see e.g. Frank, 
King & Raine 1992). In other words 



dlnar c 4 ff dlnF v i 



dlnTc dlnTc 
Using Eq. (Jl^) one obtains 



rflnT c 4 ff 
dlnT r 



din k 
dlnT r 



(18) 



(19) 
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Figure 1. Vertical structure of an accretion disc around a M = 
1OM compact object at r = 3 x 10 10 cm. a 0.1, M w 10 16 
g s" 1 (i.e. T e g « 5700 K) and T irr = 0. Both the disc temper- 
ature T and the temperature corresponding to the viscous flux 
(-Fvis/cr) 1 / 4 are plotted in units of 10 5 K. At the photosphere, 
the latter gives the effective temperature T e ff . Since T; rr = 0, the 
surface temperature T(r s ) = T e g . The dashed line is the density 
in units of 10 -7 g cm -3 and the dotted line is the ratio of the 
convective to the total fluxes (between and 1). For these pa- 
rameters, the section of the disc lays on the lower stable branch. 
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Figure 2. Same as Fig. |l|but for T irr = 6000 K. Since the irradia- 
tion temperature is non-zero, the temperature at the photosphere 
is not equal to the effective temperature (see Eq. 27). Note the 
increased disk height (same coordinate system as in Fig.[j]). 



In a gas pressure dominated disc F v [ s ~ T c (see e.g. Frank et 
al. 1992). The thermal instability is due to a rapid change 
of opacities with temperature when hydrogen begins to re- 
combine. At high temperatures dln«;/dlnT c ~ —4. In the 
instability region, the temperature exponent becomes large 
and positive dln«/dlnT c ~ 7 — 10, and in the end cooling 
is decreasing with temperature. As can be seen from Eq. 
(p"9|), even a moderate T\ rr /T c ratio modifies the transition 
from stable to unstable configurations, thereby pushingit to 
lower temperatures. This effect is clearly seen in Figs. |fj|and 
Q One should also note that another effect of irradiation 
is to suppress or alter convection, even for moderate values 
of Ti rr /T c as can be seen in Figs. 1-3 (and discussed be- 
low). This effect is moderate, as the convective flux can be 
small compared to the radiative flux; it cannot be included 
in the simple analytic model developed here for illustrative 
purposes, but is taken into account in our numerical calcu- 
lations. 



3 SELF-CONSISTENT MODELS OF X-RAY 
IRRADIATED ACCRETION DISCS 

Having demonstrated above that the simplest models taking 
into account X-ray irradiation of accretion discs in LMXBs 
give interesting results, we decided to undertake more rig- 
orous and detailed calculations to investigate these effects 
further. In all cases, we rewrite Eq. ( |l5| ) as: 



T: rr C 



with 



Mc 2 
4naR 2 



C = V (l-e) 



Hiti 



V dlnR J 



(20) 



(21) 



and the results depend on the value of C, which in turn is 
a function of r/,e and Hi T1 (R). In a self-consistent model 
Hi IT (R) results from calculations so that r\ and e are the 
only free parameters. There is not much freedom in r\ but the 
value of e is rathe r uncertain, as discussed, for example by de 
Jong et al. (1996). We assumed that X-rays are deposited at 
the disc photosphere, Hi rT = z s (R), whose height is in turn 
modified by irradiation. Our model is fully self-consistent, 
in the sense that this feedback mechanism is taken into ac- 
count. 
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Figure 3. Same as Fig. |lj but for T irr = 12000 K. Here, the 
convective flux is negligible with F conv /F to t ~ 0. The disk height 
has increased along with the irradiation temperature. 



3.1 The vertical structure 

We use the numerical code described in H98 to find the 
vertical structure of an X-ray irradiated disc. 

The equations describing the vertical structure are the 
same as in H98, with the exception of the surface bound- 
ary condition which is modified to account for illumination. 
These equations are: 

dP_ 

dz 



Tz =2p - 

dlnT d\nP, 



dz 



dz 



dF 2 3 

—7— = -a cB il K R; 
dz 2 



(22) 
(23) 
(24) 
(25) 
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where g z = Q^z is the vertical component of gravity, ? is the 
surface column density between —z and +z, a e « the effec- 
tive viscosity coefficient (see H98) and V the temperature 
gradient of the structure, which in the radiative case is 



V ra d — 



kPF z 



(26) 



4PradC# z ' 

T ra d being the radiative pressure. In the convective case, 
we use a mixing length approximation (see H98 for more 
details). 

Eq. (^i|) is the hydrostatic equilibrium equation, Eq . 
( p3|) represents mass conservation and Eqs. ( p4|) and (j25|) 
are the energy transport and conservation equations. 

The set of equations ( f^jj^ ) is integrated between the 
disc midplane and the photosphere (r s = 2/3) with the 
boundary conditions z = 0, F z — 0, T = T c , ? = at 
the disc midplane. At the disc photosphere ? = E and 



T 4 (T s )=T c 4 ff + 7t, 



(27) 



where crT e 4 ff = F v - ls . 

Figures [jl ^ and ^ present three examples of solutions 
to Eqs. (E2h2J) in which the influence of irradiation on the 
vertical disc structure is shown. One can clearly see that 
irradiation suppresses convection and flattens the temper- 
ature profile, but even when Fi TI ~ 20T v i a the disc is not 
isothermal. One should also note the growth of the mid- 
plane temperature with increasing irradiation temperature, 
from ~ 14000 K in the non-irradiated case to ~ 25000 K for 
T irr = 12000 K. 

In a geometrically thin disk model the radial structure 
can be separated from the vertical structure calculations. 
However, one needs to evaluate the heating and cooling 
fluxes for each disk annulus. In H98, T e ff and T v j s are deter- 
mined for a number of values of T c , E, and r with the inter- 
mediate values calculated by linear interpolation. In the case 
of an irradiated disc one must add T rr and the photosphere 
height z s as input and output parameters. 

Unfortunately, this is computationally intensive so that 
in order to make the problem tractable we use only 9 equally 
spaced values of Ti rr ranging between and T c (isothermal 
irradiated disk) for each {T c ,E,r}. For a given set of pa- 
rameters {T c , E,r, Ti rr }, we use a linear interpolation in T c , 
E, and r, and a cubic spline interpolation in Tj rr in order 
to determine T e g, F v i s and z s . Despite the low number of 
points in Ti rr , the fit is usually good to within a few per- 
cent. Occasionally, it can be off by 50% during the hydrogen 
ionization transition but as pointed out in H98, even in a 
non-irradiated disk, cooling fluxes are not known more ac- 
curately. 



3.2 S-curves 

It is customary to represent stationary disc solutions at a 
given radius in a E - TorE — M plane (S-curves). For a 
given {r, a, T rr } we now find a thermal equilibrium solution 
{T c , E} such that aT 4 ff = T vis (or a = Q c ff, see H98). 

Figures ^ and ^ show disc equilibria at r — 3 x 10 10 cm 
for two values of the compact object mass (1.4 and WMq) 
and 8 values of T rr ranging from to 12000 K. To illustrate 
the effect of irradiation we show four E — T diagrams rep- 
resenting the surface temperature, the midplane ( "central" ) 
temperature and the effective temperature (or mass transfer 



rate). The disappearance of the S shape with increasing T rr 
is clear. For Ti rr <; 10000 K the midplane temperature is 
always larger than ~ 10000 K so that the thermal-viscous 
instability, corresponding to the intermediate part of the S, 
is suppressed by irradiation. 

Figures ^ and ^| are obtained using a fixed, constant ir- 
radiation temperature. In actual discs, the irradiation flux 
varies with the radius and inner mass accretion rate, so that 
irradiation modifies the "radial" disc structure. Here, we are 
interested in "self-irradiation", i.e. the irradiation of the ac- 
cretion disc by X-rays produced by the accretion flow itself. 

El-Koury & Wickramasinghe (1998) have recently con- 
sidered the effect of X-ray irradiation on accretion disc S- 
curves. Our calculations correspond to their 7=1 case. 

3.3 The radial structure 

The irradiated radial disk structure is now computed using 
the time- dependent model of H98 to calculate stationary so- 
lutions in which we have zeroed the time derivatives. The 
disk evolution will be discussed in a future paper. Results 
are shown in Figures ^| and 0. We first used Eq. (^) for the 
irradiation temperature wit h valu es of albedo, efficiency and 
H 1Tr /R from de Jong et al. ( |l996 ). Results are shown in Fig- 
ures |^ and 0. We have checked that the results depend basi- 
cally on the average value C which is ~ 5 x 10~ 4 in this case. 
In other words, the small radial de pende nce in the illumina- 
tion formula used by de Jong et al. ( 1996| ) i.e. Hi IT /R cx R 2 ^ 7 
can be neglected. For this value for C, illumination keeps the 
disk on the hot branch at a significantly larger radius than 
in the unilluminated case, i.e. the stabilizing effect is strong. 
However, if one uses Hi IT = z s , or Hi lr = H to calculate C 
self-consistently, the resulting disc structure is not affected 
by irradiation, regardless of the values used for the albedo 
and efficiency : the outer disc regions that could be modified 
by irradiation are screened by the inner disc regions. The ra- 
dial profiles of midplane temperature, surface density, etc. 
are, in Figs. |^ and [?], indistinguishable from the continuous- 
line profiles of these quantities representing a non-irradiated 
disc. 

Contrary to the often used assumption of a concave disc, 
a self-consistent disc model produces a convex disc (see Fig. 
^| [?]; Tuchman et al. (199C). In the inner (radially), hot re- 
gions of a stationary disc, where the mean opacity is well 
described by the Kramers' formula (k ~ pT~ 3 5 ) the disc 
height varies as R 9 ^ 8 , but, as temperature decreases with ra- 
dius, the opacity gradually changes and this flattens the disc 
shape. Finally, at effective temperatures & 10 4 K, opacities 
decrease with temperature and d\n z s /d\n R becomes nega- 
tive and the outer disc regions are screened from X-ray irra- 
diation. Screening becomes effective when din z s /dln R — 1, 
which occurs in regions which are still stable, as can be seen 
from Figs. ^ and [j]. Of course, the same opacity behavior 
leads to the disc thermal instability, the marginally stable 
and unstable outer disc regions are therefore screened from 
X-rays (see also Cannizzo 1998). 

The abrupt change of slope of T c , E and z s /R profiles 
which is prominent in Figs. ^ and ^ results from the drastic 
change of opacities due to hydrogen recombination. In the 
non-irradiated case this happens at a midplane temperature 
T c 40000 K. Irradiation stabilizes the disc and the critical 
midplane temperature is ~ 20000 K. Note that at these 
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Figure 4. E - T surtace , S - T c , E - T e g and S - M curves for r = 3 ■ 10 10 cm, M = 1.4M , a = 0.1, and T irr = [0, 3, 6, 9, 12] X 10 3 K. 
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Figure 6. Radial profiles of the midplane temperature, surface density, and photospheric height to radius ratio for an un-irradiated 
(continuous line) and irradiated disc (dashed line) around a 1.4 Mq compact object. In the temperature diag ram the continuous line 
represents the effective temperature. The accretion rate is M = 10 18 g s _1 , a = 0.1. Tj fI is taken from Eq. (151), with r\ = 0.1, e = 0.92, 
H iTI /R = 0.2(R/R out ) 2 / 7 . Re gions beyond the radius at which a break in the T c or S curve is seen are unstable. 
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Figure 7. Same as in Fig. H but for M = 10 Mq 
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Figure 9. Same as Fig. || but for M = 10 M 



temperatures the disc is still far from being isothermal. The 
effect of irradiation is to extend the hot disc temperature 
profile (~ R~ 3 / A ) to larger radii. Note also the difference 
between the values of z s /R for 1.4 Mq and 10 Mq compact 
objects due to the stronger gravity (at a given radius) in the 
case of a larger mass central body. 

Figs. ^ and ^ showed disc structures for a given accre- 
tion rate [M = 10 18 g s _1 ). In Figs, [s] and ^ we present 
surface density profiles for 4 values of the accretion rate, for 
a M — 1.4 Mq and M = 10 Mq compact object respectively. 
As expected, the change in the radial extension of the sta- 
ble, hot branch is larger for higher mass transfer rates, i.e. 
stronger irradiation. 

The conclusion that a planar (non-warped) stationary 
disc irradiated by a point X-ray source powered by accre- 
tion is unaffected by irradiation is unavoidable. Since obser- 
vations suggest that accretion discs in LMXBs are affected 
by X-irradiation, either the assumption of a planar disc or 
that of a point source (or both) are wrong. Accretion discs 
in stationary LMXBs could be warped due to the Pringle 



irradiation instability (Pringle, 1997). In the framework of 



AGN accretion discs, a point source above the disc is of- 
ten invoked to explain the data but its physical meaning is 
unclear. 

In what follows we use Eq. (^) with C constant, which 
should be a reasonable approximation of more accurate ir- 
radiation laws as discussed above and as argued in the next 
section. 



3.4 Stability criteria 

The hot stable solutions (slope dlnE/dlni? ~ —1) end at a 
surface density E m i n which in the unirradiated case is given 
by (H98): 



E m in — 30.82 



U0 10 cm J 



ID"' cm y S Cm " 2 ™ 

where a = 0.1 and M = 1.4 Mq was assumed. Irradiation 
extends the hot branch to lower critical surface densities at 
larger radii; the amplitude of this extension increases with 
accretion rate so that the slope of the E m i n (i?) is flatter than 
in the unirradiated case. With C constant in the irradiation 
law given by Eq. ( po| ) 

z r \ -°- 3 / b \ 8 

5^) (lO^m) m 



E| 



11.4 i 



where C is scaled by its "typical value" 5 ■ 10~ . 

In the irradiated case the critical accretion rate, below 
which no steady, stable solution can exist, is given by : 

-0.4 



Kll, « 1-5 x 10 1 



— ( 

<( 



R 



C 



10 10 cm, 

0.5 

g s~ 



,5 x 10- 4 

which can be compared to the unirradiated value of H98 

0.89 



(30) 



9.5 x 10 



15 



Mi 

Mq 



R 



lQi 



g s 



(31) 



The curves corresponding to an irradiated disc also 
show a density maximum analogous to the E max appearing 
in the non-irradiated case. According to the thermal-viscous 
disc instability model, a quiescent disc in the low state must 
satisfy E < E max everywhere, and cannot be steady, for 
any reasonable value of the mass transfer rate (see e.g. Can- 
nizzo 1993; H98). This condition is unchanged by irradiation 
because, in a quiescent disc M(R) ~ R 2S ' 5 (H98), so that 
Fvis/Firt ~ R 1 ' 65 and self-irradiation is never important. 

Eq. ( ^o| ) provides the thermal-viscous stability crite- 
rion for X-ray irradiated discs. Assuming standard relations 
(King et al., 1997a) between the disc radius and orbital pa- 
rameters of the binary one gets (Phr is the binary orbital 
period in hours) 



Meat » 2.0 x 10 15 



Mi 

M© 



M 2 

Mq 



(5 x 10- 4 ) 



g s 



(32) 



whereas King et al. ( 1997a ) get 

5/6 



M!: 



3.6 x 10 1 



Mi 
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-1/6 



pi. 33 
Jlir 
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Figure 10. Stability criterion applied to neutron star and black 
hole LMXB. The values for th e different systems have been 
taken from van Paradijs (1996). Transient systems are repre- 



sented by stars and persistent systems by triangles. The criterions 
are weakly dependent on the secondary mass and, following van 
Paradijs, we took M2 = 0.4Mq. The dotted line corresponds to 
the criterion for an un- irradiated disk (Mi = 1.4Mq), the upper 
and lower dashed line are the King et al. criterions respectively 
with and without the additional H/R factor (Eq. 33). The solid 
line is the criterion deduced here (Eq. 32). The conventions are 
identical in the SXT diagram with M\ = IOMq . 



(5 X lO" 4 ) 



(33) 



King et al. assume that C should be multiplied by an addi- 
tional factor H/R ~ 0.2 to account for the non-point source 
character of a black hole. Note that Mj^ksz is inversely pro- 
portional to C because the King et al. criterion is obtained 
simply by equating T[ TI with the temperature corresponding 
to hydrogen ionization (in th is cas e Th ~ 6500 K). 

Following van Paradijs (L996), we show in Fig. [To] the 
stability criterion in the different cases of neutron star 
LMXB and SXT (black hole LMXB). As already noticed 
by van Paradijs, the non-irradiated stability criterion can- 
not explain the important population of persistent neutron 
star LMXB. For similar values of C, our calcula tions result 
in less stable disks than what King et al. ( 1997a| ) found. Our 
result is very close to that of King et al. if their additional 
H/R factor is taken into account. However, this amounts to 
comparing different values of C and is probably coincidental. 

The critical values do of course depend on C. We use 
C ~ 5 x 10~ 4 by comparison with a formula extensively used 
in the recent literature, but, as argued above, the numerical 
values of the various parameters (such as the albedo) enter- 
ing C could be different from those usually assumed. Further- 
more, we want to stress that the assumptions (isothermal 
vertical structure) used in deriving this "standard" formula 
are not justified. The results obtained through this formula 
are reasonable probably only because it is equivalent to tak- 
ing C almost constant. It is clear that the geometry and 
physics of irradiation in disks are still poorly understood and 
that much more work is needed before a satisfying model can 
be proposed. 



4 CONCLUSIONS 

We calculated self-consistent models of X-ray irradiated 
accretion discs around stellar mass compact bodies. We 
showed that irradiation of a stationary, standard, planar 
(non-warped) accretion disc by a point-like source positioned 
near the equatorial plane cannot modify the disc structure. 
The outer disc regions, which one could expect to be mod- 
ified by irradiation, cannot intercept the X-rays because of 
self-screening. Claims to the contrary must be based on the 
use of erroneous equations and/or inadequate assumptions 
about disc properties such as shape and vertical temperature 
structure. 

We used the diffusion approximation to describe radia- 
tive transfer and a simple assumption about energy deposi- 
tion by X-rays but our conclusions are general because they 
basically result from the law of energy conservation. The 
temporal evolution of unstable irradiated disks will be dis- 
cussed in a forthcoming paper (Dubus et al., 1998). 

Observations show, however, that outer regions of accre- 
tion discs in LMXBs are irradiated so they are intercepting 
X-rays. We conclude, therefore, that these discs could be 
non-planar (i.e. warped), and/or that the irradiating source 
could be not in the equatorial plane and/or could be not 
point-like (see e.g. Dumont & Collin- Souffrin (1990)). These 
types of irradiation can be summarized in a parameter for 
which we have few theoretical constraints; in view of the un- 
certainties in the models, observations will bring important 
constraints to this problem. 
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